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Abstract 

The Feller process is an one-dimensional diffusion process with linear drift and state-dependent 
diffusion coefficient vanishing at the origin. The process is positive definite and it is this property 
along with its linear character that have made Feller process a convenient candidate for the mod- 
eling of a number of phenomena ranging from single neuron firing to volatility of financial assets. 
While general properties of the process are well known since long, less known are properties related 
to level crossing such as the first-passage and the escape problems. In this work we thoroughly 
address these questions. 

PACS numbers: 89.65.Gh, 02.50.Ey, 05.40.Jc, 05.45.Tp 



jaume.masoliver@ub.edu 


josep.perello@ub.edu 





1 



I. INTRODUCTION 



Diffusion processes are Markovian random processes with continuous sample paths. From 
a mathematical point of view they are characterized, in one dimension, by two functions: 
the drift, f(y,t), and a positive defined diffusion coefficient D(y,t) > 0. The sample paths 
of any diffusion process can thus be pictured as the continuous trajectory resulting from 
the superposition of a deterministic evolution, governed by f(y,t), and fluctuations around 
it, the latter determined by D(y,t). Denoting the process by Y(t), the diffusion picture 
becomes apparent by the fact that the time evolution of Y(t) is ruled by the stochastic 
differential equation 

dY(t) = f(Y(t),t)dt+y/D(Y(t),t)dW(t). 

where W(t) is the Wiener process, that is, a Gaussian process with zero mean, unit vari- 
ance and correlation function (W (ti)W fo)) = min(ti,t 2 )- In what follows all stochastic 
differentials are interpreted in the sense of Ito. 

The Feller process is a special kind of diffusion process with linear drift and linear diffusion 
coefficient vanishing at the origin pQ. The time evolution of the process is thus governed by 

dY(t) = \-aY(t) + P]dt + ky/Y^jdW(t), (1) 

where a > 0, (3 and k > are constant parameters. 

Both Feller and Orstein-Uhlenbeck processes (diffusion processes also with linear drift 
but constant diffusion coefficient) have been widely used, with a marked prominence of 
the latter, in the modeling of countless physical phenomena. Both share a linear drift, 
f(Y) = —aY + 0, which for a > results in a restoring force that, in the absence of noise, 
makes both processes decay toward the value j5. 

However, and contrary to the Ornstein-Uhlenbek process where diffusion is constant, the 
Feller process has a state-dependent diffusion coefficient, D(Y) = k 2 Y, which for large values 
of Y enhances the effects of noise while as Y goes to zero the effect of noise vanishes. Hence, 
when the process reaches the origin, the drift drags it toward the value /3. If /3 > the 
process, starting at some positive value, cannot reach the negative region which, in turn, 
renders the process always non-negative (otherwise the noise term in Eq. (jl]) would become 
imaginary). Therefore, for the Feller process the origin is a singular boundary that the 
process cannot cross. 
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A related question is whether or not the origin is accessible, in other words, whether the 
value Y = can or cannot be attained by the process. This is a crucial question in many 
practical situations and, as we will prove later, the answer depends on the particular value 
taken by a parameter which balances the values of (3 and k 2 . The problem of classifying 
the different types of boundaries appearing in diffusion processes was thoroughly studied by 
Feller himself during the 1950's and we refer the reader to the literature for a more complete 
account on the subject |2H3]. 

Possessing a state-dependent diffusion and, most importantly, the fact that the process 
never attains negative values have made Feller process and ideal candidate for modeling a 
number of phenomena in natural and social sciences. Theoretical biology was one of the 
first places where the process was, during 1970's, seriously considered [Sj. Perhaps the most 
prominent place is within the context of neurobiology in order to model firing of single 
neurons [EH9] . The Feller neuronal model is one of the so-called stochastic integrate and fire 
models which are simple models aiming to reproduce the membrane potential fluctuations. 
Experimental progresses has lead to the possibility of fitting real data to the Feller neuronal 
model among many others [TOTTE] . 

In a different context Capocelli and Ricciardi [15J considered the possibility of modeling 
biological population with the Feller process in order to include environmental randomness 
to the classic Malthusian growth rate [TH]. The approach [TBI [T7J [TS] represents in fact 
an alternative to the Lotka-Volterra models in ecosystems and the interest in this sort of 
problems is mostly focused on the extinction -that is, on the possibility of attaining the 
singular boundary Y = 0- and also on the unrestricted growth [T7J [18] . 

Financial markets is another field where the Feller process is widely used. It was intro- 
duced in 1985 to model term structure of interest rates receiving the name of Cox, Ingersoll, 
Ross (CIR) model and successfully evaluate bond prices P3JJ [20] • The process is also be- 
ing considered to provide a random character to the volatility of a given stock. Volatility 
obeying the Feller model jointly with a log-Brownian stochastic dynamics for the asset price 
evolution configure a two-dimensional diffusion market process called Heston model j2TJ [22] 
which is a rather useful model specially for option pricing [201 EI]- 

In all of the above mentioned situations susceptible to be modeled with a Feller process, 
the first-passage time events related, among others, to level-crossing and the triggering of 
a given signal are very significant phenomena for different reasons which depend on each 
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context. This is for instance the case of the neuronal activity since a spike generation is 
due to the crossing of a threshold by the membrane potential signal. Additionally, popu- 
lation extinctions of any type or volatility bursts in financial markets are also important 
phenomena to model and study. First-passage time is, however, a difficult topic [231427] . To 
our knowledge, for the Feller process this crucial facet has been scantily studied and only 
partially solved some years ago in the context of single neuron firing [El |8] . It is our main 
objective to address the first-passage time properties of process ([!]). 

This paper is organized as follows. In Section [III we introduce the general properties of the 



unrestricted probability density of the Feller model. Section III is devoted to the derivation 
of the first-passage time and escape probabilities with special attention to a couple of specific 



situations. Section IV is mostly focused on the derivation of the mean first-passage time. 
We finally summarize the results obtained in Section |V| 



II. GENERAL PROPERTIES OF THE FELLER MODEL 

Before addressing the main issue of this paper, let us briefly review the main traits of the 
process and the role of the boundary at the origin. 

For the rest of the paper in turns out to be convenient to scale time and the process itself 
in the following way (recall we have assumed a > 0) 

t' = at, X = ^Y, (2) 

so that the Langevin equation reads 

dX(t') = ~[X(f) - 0]dt' + y/2X{V)dW{t'), (3) 

where 9 is the only free parameter left. Its relation to /3 and k is 

0=!>O. (4) 

This parameter is called "saturation level" or "normal level" and it is the value to which the 
Feller process X(t') is attracted to. As we will shortly see, 9 has a key role in the behavior 
of the process. 

Let p(x,t'\xo) be the probability density function (pdf) for process ^ to be in state x 
at time t'\ 

p(x, t'\x )dy = Prob{a; < X{t') <x + dx\X(0) = x }. 
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This density satisfies the (forward) Fokker-Planck equation (FPE) (as long as there is no 
confusion we will drop the prime in the time variable) 

t = ^-m + ^P), (5) 

with initial condition 

p(x,0\x ) = 5(x - x ). (6) 

Recall that x = is a singular boundary of the process and no "particle" can either leave 
or enter through this boundary (see Sect. [T]). A sufficient condition for this to happen is 
that the probability flux of the process through x = is zero [3]. We will thus search for 
solutions of the initial-value problem - (Jsj) that meet such a condition, that is, 

lim { (x - 9)p(x, t\x ) + ^- \xp(x, t\x )] 1=0. (7) 
z->o [ ox J 

The expression for the pdf of the process p(x,t\xo) was first obtained by Feller himself 
many years ago using a tortuous procedure which involved the solution of a rather clumsy 
integral equation pp. In the Appendix [A] we present a simpler and more direct derivation 
based on the Laplace transform of the problem ({5])-(|6]). The final expression reads 

1 /xe-*\ (e ~ 1)/2 f x + x a e- l \ (2yfxx~tfF l \ 
p{x, t\x ) = — exp <^ — } I -i — — , (8) 



x 

where Iq-\{z) is a modified Bessel function defined as 

n=0 ' 

From Eq. (|8j we easily get the stationary pdf of the process defined as 

p st (x) = lim p(x,t\x ). 



t—>oo 



Indeed, taking into account that (cf. Eq. ([9 



Io-i(z) = ^ f) (z/2) e - 1 [l + 0(z 2 )], (10) 



and from Eq. ([8]) we obtain the Gamma distribution: 
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Pst(x) = fff/~ le ~ X - (11) 



Another property that we can easily establish is the behavior of the probability distribu- 



tion at the singular boundary located at x — 0. In effect, using Eq. (10) we see from Eq. 
g that 

e -x e~ t /(l-e~ t ) 

p{xAxq) = r(6)(i - e-r ^ 1 [1 + 0{x)] ' 

from which it follows 

9 < 1 

p(0,t\x ) = { (12) 
9 > 1, 



and 



e ~x e '/(l— e ') 



p(0,t\x ) = Zf — (6 = 1). (13) 

1 — e 1 

We thus see that when 9 > 1 the probability for the Feller process to reach the value 
x = is zero but when 9 < 1 this probability is greater than zero. In other words, if 9 < 1 
the origin is an accessible boundary, while if 9 > 1 it is not [U E] . 



III. FIRST-PASSAGE AND ESCAPE PROBABILITIES 

After reviewing the main traits of the Feller process we will now focus on level-crossing 
problems -a collective name embracing questions such as hitting, first-passage, escape and 
extreme values, among others- for that process. According to whether we are dealing with 
one-sided or two-sided barrier problems, we separate level crossing into two different issues. 
In one of them, the hitting or first-passage problem, we deal with the time that the process 
reaches some "critical" value, or "threshold", for the first time. The second issue, albeit 
closely related to the first one, concerns the time when the process first leaves a given 
interval. This is the so-called escape or exit problem. 



A. First-passage probability 

Let us first address the first-passage problem for the Feller process. The problem is solved 
when one knows the first-passage probability to threshold x c . Let us denote by W c (t|x) the 
probability of first reaching x c > when the process starts at t = from the value x > 0. 

As is well known [U |23j |27] the first-passage probability satisfies the backward Fokker- 
Planck equation 

dW c n,dW c d 2 W c 

-W = - {X - 6) ^X- +X ^ (14) 
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with initial condition 



W c (0\x) = 0, 



(15) 



and boundary condition 

W c (t\x c ) = 1. (16) 



The difficulty of solving the initial-boundary problem (14)- (16) is decreased by taking 
the time Laplace transform, 



oo 

sti 



W c (s\x) = / e- st W c (t\x)dt, 
Jo 

which reduces the original problem to the solution of an ordinary differential equation (the 
Kummer equation [28]): 



with boundary condition 



W c (s\x c ) = ~ s - (18) 



Since W c is a probability it is obvious that any solution of the problem must be finite and 
non negative for all x > 0. 



The general solution of the Kummer equation (17) is 



W c (s\x) = AF(s, 6, x) + BU(s, 9, x), (19) 

where A and B are arbitrary constants and F(s,6,x) and U(s,8,x) are the confluent hy- 
pergeometric functions of first and second kind [2H] respectively defined by 



*■(.,»,«)= ^E^^ (20) 



and 



u ( s > 9 > x ) = T (3+i-e) F{s ' 9 > x) + I ^r x9 ~ lF{s + 1 " e ' 2 " e > x) - (21) 

In order to proceed further we need to specify whether the initial value x is above or 
below the threshold x r : 
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FIG. 1. First-passage probability VF c (t|a;) given by Eq. (23) as a function of the scaled time and 
for two different parameters 9 and initial values x. Inverse Laplace transform is obtained with the 
Stehfest algorithm. 

1. Initial value below threshold (x < x c ) 



In this case x can be arbitrarily small and taking into account that (see Eqs. (20)-(21)) 



r(i-Q) 
r(s+i-0) 



9 < 1 



1 oo 9 > 1 



we see that the solution to the problem staying finite for any initial position between the 
origin and x c and for any positive value of the parameter 9 is 



W c (s\x) = AF(s,9,x) 



The boundary condition (18) fixes the value of A and 



= M&r {x £ (23) 

Figure [T] shows the numerical computation of this expression in the original Laplace space. 
We have used the well-known Stehfest algorithm [29] and the results does not show any 
computational problem. As expected the closest to the threshold the fastest de first-passage 
time probability decays. And a greater 9, corresponds to a smaller W c (i|x) as well. 
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2. Initial value above threshold (x > x c ) 



In such be arbitrarily large. Hence, taking into account that 



lim F(s, 9, x) = oo, 



while U(s,9,x) stays finite for all positive values of x [2H], we see from Eq. (19) that the 
general solution of the problem which remains finite for all x > is 

W c (s\x) = BU(s,9,x), 



and from the boundary condition ( 18 ) we conclude 

m / I x U(s,9,x) 
WJs\x) - 



(x > x c ). 



(24) 



sU(s,9,x c y 

Numerical inversion of this result is again easy to compute with standard algorithms. The 



small difference lies on the fact that the hypergeometric function of second kind (21) is 



slightly more complicated than the hypergeometric function of first kind (20) 



B. Reaching the origin 

Another interesting quantity is the first-passage probability to threshold x c = 0; that is to 
say, the probability of first attaining the singular boundary of the process. This probability 
is relevant in the firing of neurons and it was addressed some years ago by by Capocelli and 
Ricciardi |6j (see also the work of Laska et al [8]) . Let us denote by Wo the first-passage 



probability to the origin. Since our process is always positive, x > 0, we must use Eq. (24) 



in order to evaluate Wo- Setting x c = in Eq. (24) and using Eq. (22) we obtain 




W (s\x) = { 31 ^ - ' 6<1 (25) 

9 > 1. 



We will now proceed to invert Eq. (25 ) thus obtaining the first-passage probability 
in real time, something that seem to be unfeasible for any value of the threshold x c , at least 
exactly (more on this below). 

Using the property [28] 

U(s, 9, x) = x x - e U(s + 1-9,2-9, x), 



we write for 9 < 1: 



Wo(s\x) 



ir(s + i 



-x 1 - e u(s + i-e,2-e,x) 



s r(i 

which, after using the following integral representation of the Kummer function U 

1 



reads 



U(a, c, x) 



W (s\x) 



T(a) 



e -xz z a-U 1 + z \c- a -l d 



X 



1-6 



sva 



Therefore, 



W (t\x) 



X 



T(l 



e -xz z -e 



Vr 1 



l + z 



dz. 



s \l + z 



dz, 



where £ 1 stands for Laplace inversion. Since [30 



c- 1 



Q(t - a) 



where G(-) is the Heaviside step function, we have 



C 



-i 



e 



t + ln 



l + z 



eu- 



l-e- 



Hence, 



or, equivalently, 



W (t\x) 



W (t\x) 



X 



r(i - e 
i 



e- xz z-°dz, 



e-*/(l-e-*) 



r l-e 



' 1-e- 



where T(a, is the incomplete gamma function (28] 



T(a, z) 



y a - x e- y dy. 



Finally, 



W Q {t\x) 



r(i- 




rfi-0 



' l-e" 



< 1 
> 1. 



(26) 



(27) 



We remark (as shown already in Eq. (25)) that when 6 > 1 the first-passage probability to 
the origin is zero in agreement with the fact, pointed out in Sect. [TT], that if 9 > 1 x = 
is unattainable. The probability of first reaching the origin is represented in Fig. [2] where 
Wo(t|x) is shown as a function of time and for two different values of parameter model 9 
and initial value x. 
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=0.2, x=0.3 



FIG. 2. First-passage probability Wo(i|x) given by Eq. (27) as a function of time (t > 1) for two 
different values of 6 {9 < 1) and the initial values x. Note that smaller values of 8 and x imply 
bigger first-passage probability. 



As t — > oo and for 9 < 1 the first-passage probability Wo becomes equal to 1, as it is 
expected since if 9 < 1 crossing the origin is a certain event as time grows. It is, however, 
interesting to see how Wq approach the unity. To this end we use the following series 
expansion of the incomplete Gamma function [28] 



T(a,z) = r(o) - J2 



n=0 



n\ a + n 



In the limit t — > oo we then get 



r i - 0, 



xe 



1-e" 



rn 



1-0 



_ x i-e e -(i-6)t + ( e -(2-e)t) 



Therefore, 



W (t\x) = 1 - 



r(2 



_ x i-e e -(i-e)t + Q ( e -(2-«)t) 



<1). 



(28) 



An interesting expression combining a power-law in x and an exponential decay in time. 



C. Large threshold 

We will now study in some detail the interesting case of a large value of the threshold 
which is the opposite case considered above. It is clear that as x c — > oo the threshold 
becomes unreachable and the first-passage probability approaches zero. Let us see how is 
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the limiting process. These results can be a really interesting for instance when one wants 
to control financial asset volatility. A large volatility corresponds to wild fluctuations in the 
asset price. 

In the case under consideration the initial position is always below threshold and the 
starting point of the analysis must be Eq. (J23j) : 

W C {S\X) = r. 

sF(s,9,x c ) 

Since now x c — > oo we use the following asymptotic expansion of the Kummer function F 

rim . r 

1 + 



F(s,9,x c ) = l^e^xT 



and 



W r (s\x) 



r(«) 

sT(6) 



r s) 



e Xc x s c - 9 F(s,9,x) 



1 + 



Our next step is the use of the following integral representation of F 

1 f°° 

F(s,9,x) = —- e- z z s - 1 F(9,xz)dz, 

r(s) J 

where F(9,xz) is the following hypergeometric series [2H] 



F(9,xz)=J2 



n=0 



1 (xz) 
(9) n n! 



(29) 



Hence, for large values of the threshold, we approximately have 

1 



WJsx) 



sT(9) 



e~ Xc x" 



e z z s F(9,xz)dz, (x c ->• oo). 



We now proceed as in Sec. HI B The Laplace inversion of the last equation reads 



W c (t\x 



~ e Xc x e r . 



Y{9) 



Since 



oo e - Z 



F(9,xz)C 



-i 



exp 



s In 



1 (- 

s \x c 



dz. 



[x, -> oo . 



then, recalling Eq. (|26|), we have 

1 / z 



-i 



Therefore, 



C 



W c (t\x) 



s \ x 



e 



t + \n 



[z - x c e *) . 



r(9) 



x n e 



-F(9,xz)dz, (x c — y oo) 



(30) 
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Using Eq. ( 29 ) we can give an alternative expression to Eq. ( 30 ) which is somewhat more 
convenient for numerical work. It reads 

OO 

W e (t\x) * ^7n\ X c e ~ XC E IffT^^ X ^)> ^ ~+ ( 31 ) 

where T(n, x c e - *) is the incomplete Gamma function. This expression is particularly suited 
for small values of the initial position. Thus, for instance, when i = 0we write 

W c {t\0) ~ ^—x e c e- x *E x (x c e-*) , (x c -> oo), (32) 



oo e - z 



where 

Ei{x) = / dz 

is the exponential integral. 



z 



D. The escape probability 

We close this section by briefly addressing the escape problem which, as mentioned before, 
is closely related with the first-passage problem studied above. The problem at hand consists 
in knowing whether or not the process X(t), starting at some point inside an interval (a, b), 
has left this interval for the first time. The answer lies in the knowledge of the survival 
probability S ab (t\x) defined as the probability that, starting at x G (a, b), the process at 
time t has not left the interval at that time or during any previous instant of time: 

Sab(t\x) = Probja < X(t') < b; < t' < t \ a < x < b), 

where x = X(0) is the starting point. The escape probability, i.e., the probability that at 
time t the process has exited the interval (a, b) for the first time, is then given by 

W ab (t\x) = 1 - S ab (t\x). 

As is well known jH EZj the survival probability obeys the backward Fokker-Planck 
equation 

dS a b , n -,dS a b d 2 S a b 
— (X — u)— h X- 



dt dx dx 2 

with initial and boundary conditions given by 

S ab (0\x) = 1, S ab (t\a) = S ab (t\b) = 0. 
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Hence, the escape probability is the solution of the initial and boundary value problem 



(compare with Eqs. ( 14 ) (pL6 



dW, 



ab 



-(x — 9) ^^ ab 4 .,• 



d 2 W, 



ab 



dt ' dx dx 2 ' 

W ab (0\x) = 0, W ab (t\a) = W ab (t\b) = 1. 



(33) 
(34) 



Following the same reasoning as before (see Eqs. (17) and (18)) we see that the time 



Laplace transform of the escape probability VF(s|x) satisfies the boundary value problem 



x - 



d 2 W ab , dW, 



(x-6)- 



ab 



dx dx 
W ab {s\a) = W ab {s\b) 



sW ab = 0, 
1 



(35) 
(36) 



Again, The general solution of the Kummer equation (35) is 



W ab {s\x) = AF(s, 9, x) + BU(s, 9, x) 



where A and B are arbitrary constants and F and U are defined in Eqs. (20) and (21) 



Boundary conditions (36) determine the value of A and B and, after routine algebra, the 



final result for the escape probability reads 

[U(s, 9, b) - U(s, 9, a)] F(s, 9, x) - [F(s, 9, b) - F(s, 9, a)] U(s, 9, x) 



Wahislx) 



(a < x < b). 



s [F(s, 9, a)U(s, 9, b) - F(s, 9, b)U(s, 9, 



(37) 



IV. LONG-TIME ASYMPTOTIC BEHAVIOR AND MEAN FIRST-PASSAGE 
TIMES 

In the previous section we have solved the hitting and escape problems for the Feller pro- 
cess by means of the evaluation of the first-passage and exit probabilities. We have obtained 
exact analytical expressions for the time-Laplace transform of these probabilities. Unfortu- 



nately exact inversion seems to be beyond reach except for the cases in Sec. Ill B -when the 



threshold is located at the origin- and in the following Sec. IV A -with approximate expres- 
sions suitable for long times. In this section we will also obtain two important magnitudes 
associated with the problem: the mean first-passage time, T c (x), and the mean escape time, 
T ab (x). 
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A. Long-time behavior of the first-passage probability 

Let t c (x) be the first-passage time for the process, starting at x, to reach some threshold 
x c for the first time. It is a random variable depending on each realization of the process. 
Formally 

t c (x) = mi{t\X(t) > x c ;X(0) = x < x c ) 
when the initial value is below threshold, and 

t c (x) = mi{t\X(t) < x c ;X(0) = x > x c } 

when the initial value is above threshold. 

We next relate the first-passage time with the hitting probability W c (t\x) defined in the 
previous section. Note that if r c (x) is the first-passage time, the hitting probability can be 
defined as 

W c (t\x) = Prob{r c (x) < t} 

which shows that W c (t\x) is the distribution function of the first-passage time. The corre- 
sponding probability density is thus defined 

f c (t\x)dt = Prob{t < t c (x) < t + dt}, 

and it is related to the distribution W c by 

- ^W- (38) 

The moments of this distribution are 

POD 

T n (x\x c )= / t n f c (t\x)dt, 
Jo 

(n — 1, 2, 3, . . . ) and the mean first-passage time (MFPT) is the first moment: 

T c (x) = Ti(x\x c ). 
Note that in terms of the Laplace transform 

poo 

f c (s\x) = / e~ st f c (t\x)dt, 
Jo 

the first-passage moments are 

,n d n f c (s\x) 



T n (x\x c ) = (-1) 



ds n 
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s=0 



which implies that, as long as T n (x\x c ) (n = 1, 2, 3, . . . ) exist, the Laplace transform of the 
first-passage time density has the following expansion in powers of s 

■iy 



fc(s\x) = ^ 



n=0 



III 



(39) 



On the other hand the Laplace transform of Eq. (38), along with the initial condition 
W c (0\x) =0, yields 



WJs\x) 



-f{s\x). 



By combining Eqs. (39) and (40) we then have 

^ :-D n 



W c (s\x) = y K —^s n - l T n (x\x. 



(40) 



(41) 



n=0 



expansion that furnishes the basis for the asymptotic analysis of the first-passage probability 
W^la;). Indeed, the so-called Tauberian theorems prove that the long-time behavior of a 
function g(t) is determined by the small s behavior of its Laplace transform g(s) [31] . For 



the case of the first-passage probability we see from Eq. (41) that the small s behavior of 
W c is 

W c (s\x) = - s - T c (x) + 0(s) = i[l - sT c (x) + 0{s 2 )}. (42) 



where T c (x) = Ti(x\x c ) is the mean first-passage time. Note that expansion (42) may also 
be written, within the same level of approximation, as 

WJslx) ' 



which by the Tauberian theorems 



s[l + sT c {x) + 0{s2)Y (43) 
implies that the long-time behavior of the first-passage 



probability ^(tjx) is given by the Laplace inversion of Eq. (43). That is, 

W c (t\x) ~ 1 - e ~ mx \ (t oo). (44) 

We have thus obtained the long-time behavior of the first-passage probability to threshold 
x c and see that the MFPT determines the long-time behavior of the first-passage probability. 
We will next evaluate this average time for the Feller process. 



B. The mean first-passage time 



In terms of the Laplace transform of first-passage probability W c (s|x) obtaining the 

(45) 



MFPT is straightforward. In effect from Eq. (42|) we see that 

'1 



TJx) 



lim 



W c (s\x) 
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Using the findings of Sec. Ill we know that the first-passage probability has different ex- 
pressions as to whether the initial value of the process x is above or below the threshold x c . 
Let us now treat these two cases including the special case x c = 0. 



1. Initial value above threshold (x> x c ) 



In this case (cf Eq. (24)) 



WJslx 



u(s,e,x) 



and 



TJx) = lim 



sU(s, 9, x c ) ' 

i u(s,e,x c ) - u(s,e,x) 



u{ s ,e,x c ) J' (46) 

The expansion in powers of s of the Kummer function U (s, 9, x) is presented in Appendix 
[B] where it is shown that 

U(s,9,x) = l + sU 1 (x) + 0(s 2 ), (47) 

where 

U x {x) = - 6) - J U(l,l + 9,z)dz, (48) 

and if)(z) = T'(z) /T(z) is the psi function. We note that the function, defined as the indefinite 
integral: 

J U(l,l + 9,z)dz 

cannot be reduced to another Kummer function [28J nor, to the best of our knowledge, to 
any other tabulated function. 



Plugging Eqs. (47)-(48) into Eq. (46) we finally obtain 



T c (x) = f U(l,l + 6,z)dz, (x>x c ). 



(49) 



2. MFPT to the origin 



In Sects. [TTJand III (see Eq. (25)) we have seen that when 9 > 1 the origin is unattainable. 
However, if 9 < 1 this singular boundary can be reached by the process. In this later case it 
is natural to ask which is the MFPT to the origin. The question has not only an academic 
interest but is relevant in mathematical biology where x = corresponds to the potential at 
which a neuron is fired [6]. Also in econophysics it is useful to know whether volatility or 
the interest rates can drop to zero and which is the average time expected to do so. 
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Since x > 0, the expression for the MFPT to the origin, denoted by Tq(x), will be given 



by Eq. (49) with x c = 0. Unfortunately setting x c = in Eq. (49) is not possible because 



the integral is singular at the lower level. 

We proceed as follows: start with the definition of the Kummer function U given in Eq. 



(21), use the integration rule [28] 



x 

F(a, c, z)dz = — F(a, c + 1, x), 
c 

and take into account the standard property of the Gamma function T(z + 1) = zT(z). We 
write 



U(l,l + 9,z)dz 



x - j F(i,i + e,z)dz-v(e-i)x 1 - 9 F(i~e,2~e,x). (50) 



Substituting into Eq. (49) and taking the limit x c — > + , we have 

T (x)= \im\-lf F(l,l + 6,z)dz 
Xc ^ + ^ J xc 

- T(0 - 1) \x x - F{\ -e,2-9,x)- x\- e F(l -6,2-6, x c 



Using the value of the Kummer function F at the origin, F(l, 1 + 6, 0) = 1 (cf Eq. (20)) 
we have 



T (x) = ~J X F(l,l + 9,z)dz 

-T(6-l) 
Hence, if 6 < 1 we get 

1 f x 

T (x) = -- / F(l,l + 6,z)dz-T(6-l)x 1 - e F(l-6,2-6,x), 
V Jo 



x l - e F(l - 6, 2 - 6, x) - lim (x 1 - 9 ) 

z c ->0+ V ' 



(51) 



and taking into account Eq. (50) we see that in this case the expression for the MFPT to 



the origin is given by Eq. (49) with x c = 0: 



T (x)= f U{l,l + 6,z)dz (0<1). 
Jo 



(52) 



However, we see from Eq. (51) that when 6 > 1 x\ 9 — > oo as x c — > + and the process 



takes an infinite average time to reach the origin 

T (x) = oo (6 > 1), (53) 
which confirms that when 6 > 1 the singular boundary x = is unattainable. 



18 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

x 



FIG. 3. The mean first-passage time as function of the initial value x for two different values of 



9 when the critical value x c equals to one. The figure jointly shows all MFPT (57) when x < x 



and the MFPT (49) when x > x c . Putting together the two mathematical expressions allow us 



to clearly observe the distinct and asymmetric behavior of the mean first-passage depending on 
whether x < x c or x > x c . 

3. Initial value below threshold (x < x c ) 



Now (see Eq. (23)) 



WJs\x 



F(s,9,x) 



Hence, 



TJx) = lim 

s->0 



sF(s,9,x c )' 
lF{s,9,x c )-F(s,9,x) 



F(s,e, Xc ) J- (54) 

In the Appendix [B] we show that the expansion of powers of s of the Kummer function 
F(s, 9, x) is 

F(s,9,x) = l + sF 1 (x) + 0(s 2 ), (55) 



where 



iq(x) = F(l,l + 9,z)dz. 



Substituting Eqs. (55)-(56) into Eq. (]54h yields 



T c (x) 



e 



F(l,l + 9,z)dz, (x<x c ). 



(56) 



(57) 



A result we obtained few years ago [32] in another context and using a different approach. 



This result that applies for x < x c jointly with previous Eq. (49) that applies for x > x 
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are presented in Fig. |3j We show there the marked asymmetric behavior of the MFPT 
depending on whether the initial value x is larger or smaller than the critical value x c 



C. The mean escape time 

We close this section by obtaining the time taken by the process X(t) starting at X(0) = x 
to first leave a given interval (a, b), where a < x < b. This is called the escape (or exit) time 
out of an interval, T a b(x), and is formally defines as 

T ab( x ) — inf {t\ a > X{t) > b;a < x < b\. 

The exit time is a random variable characterized by a distribution function, 

Prob{r a6 (x) < t\X(0) = x} = Probja > X(t) > b\X(0) = x}, 



which is precisely the escape probability W ab {t\x) discussed in Sec. Ill The moments of the 
exit time are thus defined by 



t n dW ab (t\x), 



(n = 1, 2, 3, . . . ) and the mean escape time (MET) is the first moment: 



T^>(x)=T ab (x) 



Proceeding as in Sec. IV A| we easily see that the Laplace transform of the escape prob- 



ability can be written as (cf Eq. (41)) 



n=0 



nl 



from which it follows that (cf Eq. (45)) 

T ab (x) = lim 



Wahislx) 



(58) 



Moreover, similarly to the first-passage problem discussed above, we can easily prove that 
the long-time behavior of the escape probability is also solely determined by the MET (see 



Eq. (44)): 



W ab (t\x) ~ 1 - e 



-t/T ab (x) 



(t ->■ oo). 



(59) 
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Plugging the expression for W^slx) given by Eq. (37) into Eq. (58) and taking into 
account the small s development of the Kummer functions U and F, as expressed respectively 



by Eqs. (47)-(48) and Eqs. (55)-(56) we obtain after lengthy but otherwise straightforward 



algebra the following expression of the MET: 



T ab (x) 



N ab (x) 

Dabix) 



(60) 



where 



and 



px pb 

N ab (x) = U(l,l + 6,z)dz F{l,l + 9,z)dz 

J a JO 

pb pa 

+ / U(l,l + 6,z)dz j F(l,l + 6,z)dz 

Jx JO 
rb px 

- / U(l,l + 9,z)dz / F(l,l + 9,z)dz, 

J a Jo 



D nb (x) 



( [F{i,i + e,z) + eu{i,i + e,z)]dz. 

J a 



(61) 



(62) 



V. SUMMARY AND CONCLUSIONS 



We have fully addressed the first-passage and escape problems for the Feller process. Let 
us now summarize the main results obtained. The process is an one-dimensional diffusion 
defined by a linear drift and a linear diffusion coefficient vanishing at the origin. Feller 
process has the property of being positive, a salient characteristic which has earned the 
process some popularity in modeling several phenomena, from neural activity to financial 
markets. 

Perhaps the best way to define the process is by means of a stochastic differential equation. 
In the dimensionless units defined in Eq. ^ of Sec. [Tl] this reads 

dX(t) = ~[X(t) -6}dt+ y/2X(t)dW(t), 

(we have dropped the prime in the time variable) where W(t) is the Wiener process and 
6 > is the saturation or normal level to which X(t) is attracted to as t increases. The 
origin is a singular boundary because the noise term vanishes there. In Sec. |TT] we have 
reviewed the general properties of the processes which were mostly obtained by Feller many 
years ago. One of these properties refer to the attainability of the origin in which the normal 
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level plays a crucial role. Thus if 9 < 1 the origin is an accessible boundary while if 9 > 1 
it is not. 



The bulk of the paper is developed in Sees. Ill and IV where the first-passage and escape 



properties of the Feller process are thoroughly analyzed. The first-passage problem refers 
to the crossing by the process of certain preassigned critical level or threshold x c while the 
escape problem refers to the departure of some interval (a,b). 

The first-passage properties are fully characterized by the hitting, or first-passage prob- 
ability, defined as the probability of first reaching the threshold x c by the process at time 
t or before. We denote by W c (t|x) this probability, where x is the initial position of the 
process. This probability depends on whether the process is initially below (x < x c ) or 
above (x > x c ) threshold. We have obtained exact expressions for the Laplace transform of 
the hitting probability, 

POO 

W c (s\x) = / e- st W c (t\x)dt, 
Jo 

which are summarized as 

WJslx) -- 



F(s,e,x) ^ 

sF(s,6,x c )> X — Xc > 



u{s,e,x) 



X > X r 



sU{sfi,x c )i — 

where F and U are Kummer functions. 

In general these expressions for the Laplace transform of the hitting probability cannot 
be inverted exactly in an analytical fashion and one has to resort to the numerical inversion. 
There are some instances, however, in which we have been able to obtain analytical expres- 
sions in real time. This is the case of hitting the origin which has a significant interest in 
the firing of neurons and also in the Heston volatility model of the financial analysis. We 
denote by ^(^l^) the first-passage probability to threshold x c = 0, we have shown that 



W (t\x) 



0, 9 > 1, 

where T(a, z) is the incomplete gamma function. If 9 < 1 Wo^l 37 ) ~~ * 1 as t — > oo. This the 
expected behavior since when 9 < 1 crossing the origin is a sure event as time grows. The 
way Wo approaches unity is explicitly given by the following combination of a power-law in 
the initial position and an exponential time decay 

W (t\x) = 1 - ^6)^^ + 0{e-^% (9 < 1). 
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Another instance in which we have been able to obtain an (approximate) expression for 
the first-passage probability in real time is when threshold is large. In such a case 



I v°° > 1 x n 

W c (t\x) ~ -_a;* e -*c^ r ( n>a;ce -*) > ( Xc 00 ) ) 

where T(n, x c e~*) is the incomplete Gamma function. 

The escape problem is completely characterized by the escape probability, W a b(t\x), de- 
fined as the probability of first leaving a given interval (a,b). It is complementary to the 
survival probability S a b- 

W ab (t\x) = 1 - S ab {t\x), 

where S a b(t\x) is the probability that the process has not exited (a, b) at time t or during 
any previous instant of time. Formally, 

S ab (t\x) = Prob{X(t') G (a, 6), < t' < t \X(0) = xe (a, 6) }. 

For the Feller process we have been able to obtain the exact expression for the Laplace 
transform of the escape probability which turns out to be more involved than the first- 
passage probability. It reads 

w = [Ujs, 9, b) - U(s, 9, a)] F(s, 9, x) - [F(s, 9, b) - F(s, 9, a)] U(s, 9, x) 

ab[SlX) ™ s[F(s,9,a)U(s,9,b) - F(s,9,b)U(s,9,a)] 

where a < x < b, and F and U are Kummer functions. 

We have next addressed the problem of the mean-first passage time (MFPT) and the 
mean exit time (MET). We have shown that in terms of the first-passage time moments 
T n (x\x c ) (n — 1, 2, 3, ... ) -of which the MFPT corresponds to n — 1, 7i(a;|a; c ) = T c (x)- the 
Laplace transform of the hitting probability reads 

W c (s\x) = - + V] { -^-s n - 1 T n (x\x c ). 
s ^-^ nl 

n=l 

The MFPT to the threshold x c is then given by 

1 



T c (x) = lim 



g - W c (s\x) 



From these expressions we have been able to obtain, in terms of the MFPT, the following 
long-time asymptotic expression of the hitting probability 

W c {t\x) ~ 1 - e - t/Tc(x) , (t ->■ oo). 
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For the Feller process this analysis has led to different results according to whether 
initially the system is placed below or above the threshold: 



TJx) 



(1/6) f* B F(l,l + 6,z)dz, x<x c , 
J* [7(1,1 + 0, z)dz, x>x c . 



The MFPT to reach the origin, Tq(x), has also been analyzed with the result 



T (x) 



f*U(l,l + 6,z)dz, 6<1, 



oo, 



> 1, 



which constitute an additional proof of the fact that when 9 > 1 the singular boundary 
x = is unattainable. 

The analysis of the MFPT can be exactly carried out for the MET. The resulting expres- 
sions relating the MET with the escape probability are formally the same as those relating 
the MFPT with the hitting probability as can be seen in Sec. |IVC Thus, for instance, 



T„h(x) 



lim 

s->0 



- W ab (s\x) 



and 



W ab (t\x) 



-t/T ab (x) 



(t -> oo) 



where W a b(t\x) and T^ix) are the escape probability and the MET respectively. In the case 



of the Feller process the explicit expression for the MET is given in Eqs. (60)-(62). 



Let us finally mention that the extension of the above results to the study of the extreme 
values attained by the process, such as the maximum and minimum values, as well as their 
application to financial time series -in particular the volatility- is under present research 
and we expect getting a number of results very soon. 
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Appendix A: The probability density function 



The solution to problem (|5])-([6]) is more conveniently addressed by its Laplace transform 
with respect to x: 

POO 

p(a, t\xo) = / e~ ax p(x 1 t\xo)dx. (Al) 
Jo 

Taking into account condition ([7]), the transformed problem ([5])-(|6]) reads 



dp . .dp 
p{o,0\x ) = e-° X0 . 



(A2) 
(A3) 



Equation (A2) is a linear partial differential equation of first order whose solution can be 



obtained by the method of characteristics [35]. In effect, let the function h(o) be defined by 



the characteristic of Eq. (A2), h'(a) = — [o(l + a)] 1 , that is 

1 + a 



h(o) = In 



a 



(A4) 



Then the solution of Eq. (A2), as can be rightly seen by direct substitution, is 



p(o,t\x ) = (l+a)- i>(t + h(a)), 



(A5) 



where ip(z) is an arbitrary function to be determined by the initial condition (A3), i.e., 



4>(h(a)) = (l + a) 6 e- axo , 



which implies, after inverting Eq. (|A4|) to write a in terms of h, that 

xoe~ z 



ijj( z ) = (1 - e - 2 )- y exp 



I - e- 



Substituting this into Eq. (A5) we finally obtain 

I 



p(a,t\x ) 



[l + a(l-e-*)] ( 



cxp 



axoe 



1 + a 1-e- 



Let us now proceed to the Laplace inversion of Eq. (A6). Calling 



a = 1 — e 



b = x e~ 



(A6) 



(A7) 



simple algebraic manipulations followed by a power expansion yield 



exp 



bo 



1 + ao 



e b / a exp 



I a(l + ao)J ^n\{\ 



{b/a) n 



+ ao) 
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Plugging into Eq. (A6) we get 



n=0 K ' 

Let us denote by £~ 1 {f(cr)} = f(x) the operation of Laplace inverting f(o) and recall 
the standard property 

C- 1 {f(aa+l)} = -e- x / a f(x/a) 
a 

and also 

( 1 ^ „n+9— 1 

c- 1 11 



CT n+e j r(n + 0)" 



Then the Laplace inversion of Eq. (|A8|) yields 

n\T(n 



n=0 

which after simple manipulations reads 



1 f p\ -(x+b)/a (vbx/a 



\2n+9-l 



n\T(n + 6) 

We recognize the series appearing in this equation as the expression of a modified Bessel 
Function. Indeed 

[z/2) 2n+ » 



iu{z) = 



where l v [z) is the modified Bessel function of order v [2H]- Therefore, 

1 /I\(^l)/2 r^su, ( 2\fbx 



which after reverting to the original notation (cf. Eq. (A7)) reads 



, . . 1 ( xe 1 \^ f x + x e *1 T ( 2\/xx e~ t 
p{x, t\x ) = — ( — — ) exp <j ) J e _! 



1 — e * V xo J I 1 — e * J V 1 — e 



which is Eq. m 



Appendix B: Expansions for F(s,6,x) and U(s,9,x) 

In terms of the Pochhammer's symbol (a) n = T(a + n)/T(a) Kummer function F is 
defined as the series [2H] 

0)n X n 



F(s,9,x) = J2 



„ (f)n n\ 
n=0 v ' 
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Since (s) = 1 and 



(s) n = s(s + l)(s + 2) • • • (s + n - 1) = s(n - 1)! + 0(s 2 ) 



we have 



1 x r ' 



n=l 

In the sum of the right hand side we make the replacement n — > n + 1 and take into account 
that (9) n +i = 0(9 + l) n , we thus write 



x ^ l x n _ l >^ 



1 X 



n=l v ' n=0 

We can easily see that (l) n = n\, hence 



n+l -y 00 ^ | rx 



1 X n 1 







E 

n=0 



(l)n Z* 



dz. 



which, after recalling the definition of the confluent hypergeometric function F, Eq. (20) 
yields 

°° 1 T n i rx 

Substituting into Eq. (|B1|) we get 



F(s, 0, x) = 1 + sFx(x) + 0(s 2 ), 



(B2) 



where 



F x (x) = \ f F(l,l + 9,z)dz. (B3) 
Jo 

The small s expansion of the Kummer function of second kind U(s, 8, x) is a bit more 
involved. We start from the definition of U in terms of F (cf Eq. (|21~j)) 



U{s,6,x) 



r(i 



F(s, 0, x) + .^ x^ 1 ^! - + s, 2 - 0, x) 



r(l -0 + sf 1 T(s 

then expand 

r(i - e + s) = r(i - 0) + sr'(i - 0) + o(s 2 ) 
= r(i-0) [i + #(i-0) + o(s 2 ) 

where ip(z) = T'(z)/T(z) is the psi function [2H]- Also [31] 

ns) = -[i- ls +o(s 2 )}, 

s L J 
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(B4) 



(B5) 



(B6) 



where 7 = 0.5772 • • • is Euler's constant. Plugging Eqs. (B2), (B5) and (B6) into Eq. (B4) 

we get 



U{s,9,x) = l + sU 1 (x) + 0(s 2 ), 



where 



Ut(x) = F 1 (x) - V(l -B) + T(6 - 1)F(1 -9,2-9,x). 
Let us finally show that a more convenient form for U\{x) is given by 

Ui(x) = - 9) - [ 17(1, 1 + 9, z)dz.. 



(B7) 
(B8) 

(B9) 



In effect, recalling the definition of F\(x) given in Eq. (B3) and using the integration rule 

f x c ~ 1 F(a, c, x)dx = (x c /c)F(a, c + 1, x), we have 

~ n- 



c/i^) = -v(i - 



|r(i 



i + e,x) + v(e)x e F{\ -e,i-e,x) 



dx. 



where we have used the well known property T(l+z) = zT(z) to write 1/8 = —T(—9)/T(l—9) 
and (1 — 9)T{9 — 1) = —Y{9). Note that the integrand is precisely the Kummer function of 
second kind 



u(i,i + e,x) 



r(- 



r(i 



F(l, l + e,x) + T(9)x e F(l -9,l-9,x) 



see Eq. (B4) with s = 1 and 9 replaced by 1 + 9). We have thus proved Eq. (B9). 
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